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general  belief  that  an  explanation  can  he  given  by  means  of  a  i>onlinear  large  deflectio: 
theory.  The  current  theories  lean  heavily  upon  a  plausible  but  arbitrarily  chosen  ener¬ 
gy  criterion  first  given  by  Tsien,  and  predict  that  the  initial  buckling  load  should  be 
influenced  by  the  testing  machine  stiffness.  Two  experiments  described,  statistically 
analysed,  and  discussed  by  the  authors  in  a  previous  paper  revealed  that  test  machine 
extensional  rigidity  does  not  influence  the  initial  failing  load.  These  results  were 
obtained  from  repeated  tests  on  a  single  near-perfect  aluminum  specimen  and  individual 
tests  on  many  less  perfect  steel  shells.  They  had  R/t  ratios  of  313  and  226,  respect¬ 
ively,  and  buckled  at  77  and  45  percent  of  the  classical  critical  load. 

Two  additional  experiments  are  described  in  detail  in  this  report  supplement. 
They  extend  the  authors’  previous  research.  Circular  cylinders  with  higher  R/t  ratios 
which  buckled  at  a  lower  percentage  of  the  classical  critical  load  were  used.  The  shellt 
had  R/t  ratios  of  946  and  1419  and  buckled  at  43.9  and  24,6  percent  of  classical,  res¬ 
pectively.  In  all  four  experiments,  204  tests  were  conducted.  Ranges  in  test  machine 
extensional  stiffness  from  589,000  lb/in. to  2400  lb/in. were  considered;  R/t  values 
ranged  from  226-1419  and  P^^/P  ,  values  from  .249  to  .77*  L/D  ratios  were  approximately 
constant ‘in  all  caseb.  The  results  show  that  the  Tsien  criterion  is  inapplicable.  This 
may  he  due  to  the  invalidity  of  the  criterion  itself  or  to  the  inadequacy  of  the  large 
displacement  analysis.  Thus,  scatter  in  experimental  results  is  not  due  to  the  influence 
of  test  machine  rigidity,  and  a  new  look  at  the  large  displacement,  analysis  and  the 
appropriate  criterion  is  reeded. 
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ABSTRACT 


The  inadequacy  of  the  classical  small  displacement  theory  of  thin  shells  in 
explaining  the  buckling  phenomenon  for  circular  cylindrical  shells  under 
axial  compression  and  spherical  shells  under  external  pressure  loading  has 
long  been  established.  Likewise,  it  has  become  general  belief  that  an 
explanation  of  their  behavior  can  be  given  by  means  of  a  nonlinear  large 
deflection  theory.  The  current  interpretations  of  such  theories  when  appl¬ 
ied  to  these  problems  lean  heavily  upon  a  plausible  but  arbitrarily  chosen 
energy  criterion  first  given  by  Tsien,  and  predict  that  the  initial  buck¬ 
ling  load  of  thin  shells  should  be  influenced  by  the  testing  uiachine  stiff¬ 
ness.  Two  experiments  described,  statistically  analyzed,  and  discussed  by 
the  authors  in  a  previous  report  revealed  that  test  machine  extensional 
rigidity  does  not  influence  the  initial  failing  load  of  axially  compressed 
cylinders  to  a  high  degree  of  probability.  These  results  were  obtained 
from  repeated  tests  on  a  single  near-perfect  aluminum  specimen  and  individ¬ 
ual  tests  on  many  less-perfect  steel  shells .  They  had  R/t  ratios  of  313 
and  226,  respectively,  and  buckled  at  77  and  45  percent  of  the  classical 
critical  load. 

Two  additional  experiments  using  the  single  specimen  approach  described 
in  detail  in  this  report  supplement  and  extend  the  authors'  previous  res¬ 
earch  by  considering  circular  cylinders  with  higher  R/t  ratios  which 
buckled  at  a  lower  percentage  of  the  classical  critical  load.  The  shells 
had  R/t  ratios  of  946  and  1419  and  budded  at  43-9  and  24.6  percent  of 
classical,  respectively.  In  all  4  experiments,  204  tests  were  conducted. 
Ranges  in  test  machine  extensional  stiffness  from  589,000  lb/in.  to  2400  lb  / 
in.  were  considered;  R/t  values  ranged  from  226  to  l4l9,  and  P  /P  values 
ranged  from  .249  to  .77*  L/*D  ratios  were  approximately  constESt  in  all 
cases.  The  results  provide  overwhelming  evidence  that  the  Tsien  criterion 
is  inapplicable  in  all  problems  considered.  This  may  be  due  to  the  invalid¬ 
ity  of  the  criterion  ■’tself  or  to  the  inadequacy  of  the  large  displacement 
analysis.  The  conse  quences,  however,  are  the  same;  scatter  in  experimental 
results  is  not  due  to  the  influence  of  test  machine  rigidity,  and  a  new  look 
at  the  large  displacement  analysis  and  the  appropriate  criterion  is  needed. 
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INTRODUCTION 


The  veracity  of  the  conmonly  accepted  Tsien  criterion  for  the  instability 
of  shell  bodies  has  been  examined  by  the  authors  in  recent  publications  .* 
According  to  this  belief— constant  total  potential  energy  before  and 
after  buckling — there  should  be  a  difference  in  the  initial  axial  buckling 
load  of  a  circular  cylindrical  shell  between  tests  made  in  rigid  and  dead¬ 
weight  testing  machines.  These  two  extreme  cases  are  illustrated  by  the 
dotted  lines  in  the  equilibrium  curve  of  the  Karman-Tsien5  post.buckling 
theory  shown  in  Figure  1.  Normal  elastic  testing  machines  are  typified  by 
the  solid  line.  It  was  pointed  out  by  these  authors  that  the  elastic 
characteristic  of  the  testing  machine  might  be  a  cause  of  the  large  scatter 
of  the  data  obtained  by  different  experimenters. 


PRIOR  RESEARCH  ON  THE  PROBLEM 


Prior  to  the  research  cited  above,  little  experimental  work  had  been 
carried  out  to  check  this  premise.  The  first  experiments  with  a  direct 
bearing  on  the  issue  appear  to  be  those  of  Horton,  Johnson,  ~nd  Hoff. 

These  tests  showed  that  the  effect  was  questionable  and  indicated  strongly 
the  need  for  a  more  intensive  program.  A  later  Btudy  by  Mossakovskii  and 
Smelyi?  appeared  to  verify  the  Karman-Tsien  theory,  but  on  close  examina¬ 
tion  the  finality  of  the  conclusion  is  marred  by  the  paucity  of  data  on 
which  it  was  founded.  Almroth,  Holmes  and  Brush®  noted  that  the  test  vehi¬ 
cle  characteristics  appeared  to  be  more  important  than  those  of  the  test 
machine;  but  this  sound  comment  can  only  be  considered  as  qualitative, 
since  no  quantitive  data  were  presented.  Likewise,  the  observation  of 
Krenzke^  that  his  repeated  tests  on  a  single  plastic  sphere  gave  no  evi¬ 
dence  to  support  Tsien' s  criterion  can  only  be  regarded  as  a  Justification 
for  further  study. 

EXPERIMENTAL  APPROACH  FOR  FURTHER  RESEARCH 


The  extensive  experimental  studies  reported  in  Reference  2  have  revealed, 
after  a  careful  statistical  reduction  of  the  data,  that  test  machine 
extensional  rigidity  does  not  influence  the  initial  failing  load  of  axially 
compressed  cylinders  to  a  high  degree  of  probability.  These  results  were 
obtained  from  repeated  tests  on  a  single  near-perfect  aluminum  specimen 
(test  series  B)  and  individual  tests  on  many  less-perfect  steel  shells 
(test  series  A).  They  had  R/t  ratios  of  313  and  226,  respectively,  and 
buckled  at  77  and  45  percent  of  the  classical  critical  load. 

It  could  be  argued  that  in  these  two  tests  the  ranges  of  R/t  and  PcrAcl 
considered  exclude  those  of  the  greatest  practical  interest.  Likewise, 
it  might  be  asserted  that  the  audio-visual  method  of  buckle  determination 
used  there  was  inadequate.  To  counter  these  possible  criticisms,  two 
additional  experiments  (test  series  C  and  D)  were  conducted  using  the 
single  specimen  approach;  they  are  described  in  detail  in  this  report. 

The  shell  specimens  had  R/t  ratios  of  9^6  and  1419  and  buckled  at  43.9 
and  24.6  percent  of  classical,  respectively. 
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Figure  1.  Equilibrium  Curve  for  Buckling  of  a  Thin  Circular 
Cylindrical  Shell  Unver  Uniform  Axial  Compression. 
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DETAIL  OF  THE  TEST  SPECIMENS 


Both  specimens  were  made  from  12-inch-wide  precision-rolled  shim  steel 
material.  One  test  vehicle  was  made  from  2-mil  stock  and  the  other  from 
3-mil  stock.  The  actual  shells  were  manufactured  by  a  wraparound  and 
seam  technique,  using  an  aluminum  mandrel  which  had  been  accurately  mach¬ 
ined  between  centers  to  a  diameter  of  5 *677  inches  with  less  them  3/10,000- 
inch  taper  in  the  full  12-inch  length.  The  flat  sheets,  cut  to  proper 
size,  we  .  lap-jointed  and  soft  soldered.  Removal  from  the  manufacturing 
mandrel  was  readily  accomplished  by  slightly  tapping  or  cooling  the  mandrel- 
specimen  assembly.  The  setup  is  shown  in  Figure  2. 

Constant  edge  restraint  was  assured  throughout  the  test  sequence  by  casting 
the  edges  of  the  cylinder  into  stiff  end  plates  with  a  low-melting  temper¬ 
ature  alloy  (Cerrolow)  which  contracts  slightly  upon  cooling.  Inward  buckle 
motion  was  limited  to  one  wall  thickness  by  means  of  an  accurately  machined 
interior  mandrel  which  was  pin  positioned  on  the  bottom  end  plate  to  remain 
concentric  with  the  shell.  The  individual  components  and  the  assembled 
unit  are  shown  in  Figure  3« 


TEST  ENVIRONMENT 


The  Baldwin-Lima-Hamilton  Universal  testing  machine  used  in  test  series  A 
was  again  employed  here.  Its  extensional  stiffness  variations  were  obtain¬ 
ed  using  the  previously  adopted  principle  as  shown  in  Figure  4(a).  It 
consisted  of  a  simple  steel  leaf  spring  system  with  a  supporting  and  clamp¬ 
ing  bracket  which  allowed  the  number  of  leaf  springs  and  their  support 
positions  to  be  varied. 

The  bracket  was  securely  mounted  to  the  upper  platen  of  the  testing 
machine,  as  shown  in  Figure  4(b).  This  permitted  the  number  of  leaf  springs 
to  be  varied  without  any  alteration  whatsoever  in  the  test  setup  of  the 
specimen  below.  In  addition,  the  lateral  and  torsional  rigidity  of  the 
basic  machine  was  not  affected.  A  magnetized  pad  of  ground  tool  steel 
positioned  on  the  lower  spring  provided  a  hard,  smooth  contact  point  for 
the  steel  loading  ball  on  the  top  of  the  specimen  which  was  identical  from 
test  to  test. 

The  stiffnesses  of  the  basic  machine  and  seven  other  leaf  spring  modifi¬ 
cations  were  determined  by  forcing  apart  the  lower  loading  platen  and 
contact  point  on  the  upper  leaf  springs  with  a  60-ton-capacity  hydraulic 
jack.  Details  of  a  typical  test  setup  are  clearly  depicted  in  Figure  5* 
Figure  5(c)  Bhows  the  quarter- span  supports  which  were  used  with  the  6  leaf 
spring  arrangement  to  provide  increased  rigidity.  In  Figure  5(d),  a  rigid 
steel  block  clamped  between  the  single  leaf  spring  and  the  supporting 
bracket  developed  the  full  rigidity  of  the  basic  machine.  Care  was  taken 
to  assure  that  the  cross-beam  of  the  testing  machine  was  positioned  at  the 
height  used  in  the  testing  sequence.  The  motion  which  resulted  and  the 
load  which  was  induced  were  measured  by  two  dial  gages  and  the  machine  load 
scale,  respectively.  In  all  cases,  the  dial  gages  were  placed  symmetrically 
about  the  load  point  in  the  plane  of  the  fixed  cross-beam. 


(a).  View  of  Mandrel-Specimen  Assembly 

Showing  Detail  of  Soft  Solder  Joint. 


(b).  Finished  Specimen  Removed  From  Manufac¬ 
turing  Mandrel. 


Figure  2.  Method  of  Making  Shim-Steel  Specimens  -  Test  Series  C  and  D. 
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(c).  View  of  Quarter-Span  Supports  Used  to 
Increase  Rigidity. 


(d).  Steel  Block  in  Position  to  Develop  Full 
Sniff ness  of  Basic  Machine. 

.  Continued . 
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The  force-head  separation  data  are  recorded  in  Tables  I  through  VIII. 

The  composite  stiffness  of  each  configuration  is  determined  from  the  load- 
average  deflection  plots  shown  in  Figures  13  through  20  and  is  summarized 
in  Table  X.  The  table  shows  stiffness  values  ranging  from  577,000  for 
7,800  lb/inch — a  74:1  variation. 

METHOD  OF  LOAD  APPLICATION  AND  DETERMINATION  OF  INSTABILITY 


Since  all  stiffness  variations  were  accomplished  by  means  of  an  overhead 
leaf  spring  system,  no  movement  of  the  test  vehicle  was  necessary  during 
the  entire  test  sequence.  However,  to  assure  uniformity  of  load  distri¬ 
bution  after  each  stiffness  alteration,  the  end  shortening  of  the  cylinder 
was  monitored  at  three  equally  spaced  points  around  the  periphery  using 
sensitive  strain  gage  deflectometers.  The  analog  signals  received 
from  all  three  transducers  were  recorded  simultaneously  on  a  Sanborn 
recorder.  A  well-lubricated  steel  loading  ban  located  on  the  top 
end  plate  was  positioned  in  such  a  manner  that  changes  in  end  shortening 
shown  by  these  readings  were  equal  to  within  2/10,000  inch  in  the  initial 
stages  of  loading. 

Buckle  determination  was  made  by  using  an  electro-optical  noncontacting 
displacement  probe  to  monitor  displacement  at  a  point  on  the  shell  wall. 
This  device,  a  Fotonic  Sensor,  uses  a  fiber  optic  cable  to  direct  a  con¬ 
stant  intensity  light  source  on  a  moving  surface,  and  it  can  detect  var¬ 
iations  in  the  amount  of  reflected  light  (Figure  6 ) .  Such  equipment  has 
good  sensitivity  and  high  resolution  and  can  be  used  to  measure  displace¬ 
ments  on  the  order  of  microns.  In  the  test  series  it  was  positioned  in 
the  following  manner.  The  shell  was  buckled  In  the  base  machine  and  shown 
to  be  invariant  in  load-carrying  capability  and  buckle  location.  The  probe 
was  then  positioned  normal  to  the  shell  wall  near  the  spot  where  initial 
buckling  was  known  to  occur.  This  procedure  is  shown  in  Figure  7.  It  wab 
subsequently  calibrated  in  place,  and  the  output  signal  was  monitored  on 
the  Sanborn  recorder  shown  in  Figure  8 . 

Figure  9  depicts  the  typical  wall  motion  which  resulted  during  loading 
of  the  3-rail  specimen.  Load  values  were  read  from  the  Baldwin-Lima- 
Hamilton  load  cell  and  were  instantly  recorded  on  the  strip  chart  by  an 
electric  impulse  marker.  The  point  of  buckling  is  clearly  defined.  A 
similar  load-radial  deflection  history  is  given  in  Figure  ID  for  the  2-mil 
specimens.  This  loading  process  was  repeated  six  times  at  each  of  eight 
different  levels  of  machine  stiffness  taken  in  a  random  fashion  for  both 
the  3-rail  and  the  2 -mil  test  vehicles.  Although  the  loading  rate  did  not 
vary  over  a  wide  range  during  the  test  series,  no  effort  was  made  to  insure 
uniformity . 


EXPERIMENTAL  DATA 

All  critical  load  values  obtained  for  the  3-mil  specimen  are  given  in 
Table  X.  Likewise,  the  buckling  loads  for  the  2-mil  cylinder  are  recorded 
in  Table  XII.  Each  test  was  conducted  according  to  the  procedure  outlined 
la  the  preceding  section. 
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The  mean  value  of  initial  failing  load  from  all  48  tests  of  the  3“®il  f^st 
vehicle  was  456.72  lb,  giving  a  ratio  of  (Fcr/Pcl)  =  0.439.  The  corre¬ 
sponding  value  for  the  2-mil  shell  was  112.59  l"b,  yielding  a  ratio  of 
(Fcr/Pcl)  =  0.246. 

DISCUSSION  OF  EXPERIMENTAL  RESULTS 


Since  several  levels  of  machine  rigidity  were  considered  in  these  last  two 
test  series,  a  regression  analysis  (adopted  previously  in  test  series  B2) 
was  again  employed  to  examine  the  possibility  of  a  relationship  between 
the  stability  level. and  the  machine  extensions!  rigidity.  Here,  an  even 
more  powerful  analysis  can  be  made  because  of  the  greater  number  of  degrees 
of  freedom.  The  lengthy  calculations  are  performed  in  Appendixes  II  and 
III,  but  the  important  results  are  presented  below. 

It  is  assumed  that  the  buckling  load  may  be  expressed  as  a  linear  function 
of  machine  stiffness.  The  method  of  least  squares  leads  to  the  empirical 
regression  lines  drawn  in  Figures  11  and  12  for  the  3-mil  and  2-mil  buck¬ 
ling  data,  respectively.  These  lines  are  drawn  through  the  discrete 
average  buckling  loads  obtained  for  each  set.  The  theoretical  slopes  {3  and 
the  intercepts  A'  are  examined  statistically  as  before  to  check  the  initial 
assumption  that  buckling  load  is  related  to  machine  rigidity.  The  same 
significance  tests  used  in  test  series  B  apply  here,  but  a  value  of  k  =  6 
must  be  used . 


In  addition,  for  these  two  cases  it  is  desirable  to  compare  the  variation 
about  the  regression  line  with  that  existing  within  the  sets  of  load  values 
corresponding  to  the  individual  machine  stiffness.  This  is  done  with  the 
F  test,  and  for  both  cases  the  hypothesis  of  linearity  is  accepted  at  the 
95-percent  confidence  level.  Additionally,  the  homogeneity  of  the  several 
set  variances  is  established  by  Cochran's  test,  and  the  test  shows  for  both 
shells  that  the  variances  do  not  differ  significantly.  Thus,  it  is  per¬ 
missible  to  make  a  pool  estimate  of  variance  about  the  regression  line  and 
to  establish  a  pooled  standard  deviation.  The  values  of  these  quantities 
are  5*3939  It.2  and  2.3225  lb.  for  the  3-mil  shell  and  4.0152  lb.2  and 
2.004  lb.  for  the  2-mil  shell.  Thus,  the  standard  deviation  of  slope  and 
intercept  can  be  computed.  These  values  are  0.1862  x  lO^in.  and  0.39*+  lb 
for  the  3-mil  shell  and  O.I606  x  10"^ in.  and  0.340  lb.  for  the  2-mil  shell. 
The  95-percent  confidence  interval  for  (3  is  [-0.5958;  0.1546]  x  10"5in. 
and  for  the  intercept  is  [455 *5^;  458.38]  lb  for  the  3-nil  specimen  and 
[-O.559,  0.088]  x  10  'in.  and  [112.17;  113*53]  It,  respectively,  for  the 
2-mil  shell.  The  t  test  statistics  to  check  the  hypotheses  of  zero  slope 
and  intercept  equal  to  the  average  critical  load  are  -I.185  and  0.623, 
respectively,  for  the  3-®il  specimen  and  -1.468  and  0.772,  respectively, 
for  the  2-mil  shell.  The  criterion  for  rejection  of  these  hypotheses  is 


I  I  *  V2i  ka  -  2 


where  k  is  the  number  of  tests  within  each  set  and  n  is  the  number  of  sets. 
At  th_  5-percent  level  of  significance  with  k  =  6  and  n  =  8, 
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Figure  12.  Regression  Analysis  for  Test  Series  D;  Buckling  Load  Considered  as  a  Linear  Function  of 
Machine  Stiffness. 


Since 


|  0.623  |  <  |  0.772  |  <  |  -1.185  |  <  j  -] .468  |  <  2.015 
all  hypotheses  are  accepted  at  the  95 -percent  level. 

The  extreme  sensitivity  of  the  foregoing  analysis  is  readily  demonstrated. 
The  operating  characteristic  curve  for  the  t  test  at  the  5-percent  level 
of  significance  and  a  sample  size  of  48  shovs  a  95 -percent  probability 
that  the  hypotheses  (the  true  slopes  are  zero)  would  be  rejected  if  they 
differ  from  zero  by  as  little  as  O.7659  x  lO"?  in.  and  O.6606  x  10"5  in. 
for  the  3-mil  and  2-mil  cylinders,  respectively.  These  slopes  imply  changes 
in  load  value  at  zero  stiffness  of  only  0.23  and  0.66  percent  of  the  corre¬ 
sponding  average  critical  load  levels.  Likewise,  the  same  curve  shows  a 
95-percent  probability  of  rejecting  the  hypotheses  that  the  true  intercepts 
are  equal  to  the  mean  values  of  buckling  load  when  they  differ  from  these 
values  by  as  little  as  0.35  and  1.24  percent  of  the  corresponding  averages. 

CONCLUSIONS 


The  two  experimental  studies  described  in  Reference  2  (teBt  series  A  and  B) 
and  the  two  presented  here  (test  series  C  and  D)  have  produced  results 
which  are  in  absolute  accord  with  each  other.  It  is  clear  from  the  several 
statistical  analyses  that  the  data  are  of  high  quality,  and  the  probability 
of  error  in  the  result  is  slight.  The  work  demonstrates  that  test  machine 
rigidity  has  no  influence  on  initial  instability  load  for  circular  cylindri¬ 
cal  shells  in  the  range  of  practical  interest.  Likewise,  from  the  work  of 
Reference  4,  system  stiffness  used  for  external  pressure  loading  of  such 
shell  bodies  does  not  affect  their  initial  crippling  pressure.  As  a  result, 
it  must  be  considered  that  the  Tsien  criterion  of  instability  is  improper 
or  that  the  large  displacement  analysis  is  inapplicable,  or  both. 

3 

These  findings,  reported  in  a  compendium  elsewhere  ,  are  in  agreement  with 
the  majority  of  earlier  predictions  based  upon  scanty  experimental  evidence 
and  qualitative  argument  but  are  antithetic  to  the  conclusions  reached  by 
Mossakovskii  and  Smelyi.  This  direct  contradiction,  however,  merely  serves 
to  emphasize  the  danger  in  statistical  deductions  made  from  an  inadequate 
volume  of  test  data. 
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APPENDIX  I* 

DETERMINATION  OF  TEST  MACHINE  EXTENSIONS  RIGIDITif -SERIES  C  TESTS 


The  rigidity  of  the  test  machines  used  in  this  investigation  was  obtained 
in  accordance  with  the  procedure  outlined  on  page  3  of  this  report.  In 
this  appendix,  the  actual  test  data  obtained  are  presented  and  analyzed. 
The  load  deflection  histories  are  presented  in  Tables  I  through  VIII,  and 
are  displayed  graphically  in  Figures  13  through  20.  The  actual  machine 
stiffnesses  are  given  on  the  appropriate  figures. 


*  See  appendixes  1  and  2  of  Reference  2  for  additional  information. 
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TABLE  I.  FORCE-HEAD  SEPARATION  DATA  FOR  THE  60,  OOO-LB  BALDWIN-LIMA- 
HAMUTON  TESTING  MACHINE 


Load 

(lb  ) 

Deflection,  Deflection, 

Dial  No.  1  .  Dial  No.  2  . 

(inches  x  10  )  (inches  x  10  ) 

Average 

Deflection  . 

(inches  x  10  ) 

85 

1-5 

1 

1.25 

140 

2.5 

2 

2.25 

195 

3.5 

3 

3.25 

265 

4.5 

4 

4.25 

350 

6.5 

6 

6.25 

385 

7 

6.5 

6.75 

475 

9 

8 

8.5 

535 

10 

9 

9.5 

610 

12 

10 

11 

755 

13.5 

12 

12.75 

860 

15.5 

14 

14.75 

1145 

20.5 

19 

19.75 

1240 

23 

20.5 

21.75 

1440 

26 

23 

24.5 

Plot  of  data  Is 

shown  in  Figure  13. 

_ ! 
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STIFFNESS  CONFIGURATION 
BASIC  MACHINE 

BASIC  STIFFNESS 


Figure  13 •  Stiffness  Plot  of  the  Basic  60, 000- lb  Baldwin-Lima' 
Hamilton  Testing  Machine  -  Test  Serins  C  and  D. 
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TABLE  II. 

FORCE-HEAD  SEPARATION  DATA  FOR  THE  60,  OOO-LB  BALDWIN-LIMA- 
HAMUTON  TESTING  MACHINE  MODIFIED  BY  6  LEAF  SPRINGS  WITH 
QUARTER-SPAN  SUPPORTS 

Load 

(lb  ) 

Deflection, 

Dial  No.  1  ^ 
(inches  x  10  j 

Deflection, 
Dial.  No.  2  . 
(inches  x  10  j 

Average 

Deflection  . 

(inches  x  10  ; 

55 

10 

0 

5 

110 

17.5 

0 

8.75 

225 

?4 

6 

15 

280 

29 

9 

19 

390 

36.5 

17 

26.75 

470 

41.5 

24 

32.75 

550 

47.5 

31.5 

39.5 

610 

50 

37 

43.5 

700 

56.5 

46.5 

51.5 

760 

59-5 

52 

55-75 

870 

67 

62.5 

64.5 

960 

71 

71 

71 

1070 

77 

8O.5 

78.75 

1200 

82 

92 

87 

I3to 

89 

104.5 

96.75 

1490 

97 

115.5 

106.25 

'lot  of  data  is  shown  in  Figure  34 


LOAD  *10" 


COMPOSITE  STIFFNESS 
=138,000  Ib/inch 


Figure  1*)-. 


Stiffness  Plot  of  bO, 000-lb  Baldwin-Lima-Hamilton 
Testing  Machine  Modified  by  Leaf  Spring  Configur¬ 
ation  62  (With  Quarter-Span  Supports)-  Test  Series 
C  and  D. 
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TABLE  III. 

FORCE-HEAD  SEPARATION  DATA  FOR  THE  60, OOO-LB  BALDWIN-LIMA- 
HAMILTON  TESTING  MACHINE  MODIFIED  BY  6  LEAF  SPRINGS 

L  1 

Load 

Deflection, 

Deflection, 

Average 

Dial  No.  1  . 

Dial  No.  2  . 

Deflection  . 

(lb  ) 

(inches  x  104) 

(inches  x  10  j 

(inches  x  10  ; 

6o 

14 

8 

11 

130 

27.5 

19 

23.25 

185 

38 

32 

35 

245 

46 

44 

45 

300 

58 

58 

58 

360 

67 

73 

70 

450 

84 

95 

89.5 

555 

101 

115 

108 

635 

117 

131 

124 

800 

148 

162 

155 

920 

170.5 

186 

178.25 

1150 

213 

232 

222.5 

1200 

224 

244 

234 

1290 

240.5 

260.5 

250.5 

1450 

272 

290.5 

281.25 

E 


'lot  of  data  is  shown  in  Figure  15 


COMPOSITE  STIFFNESS 
=  51,200  Ib/inch 


Figure  15.  Stiffness  Plot  of  60,000-'lb  Baldwin-Lim-ilamilton 

Testing  Machine  Modified  by  Leaf  Spring  Configuration 
ol  -  Test  Series  C  ana  D. 
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TABLE  IV. 

FORCE-HEAD  SEPARATION  DATA  FOR  THE  60,000-LE  BALDWIN-LIMA- 
HAMILTON  TESTING  MACHINE  MODIFIED  BY  5  LEAF  SPRINGS 

Load 

(lb  ) 

Deflection, 

Dial  No.  1  j, 
(inches  x  10  ) 

Deflection, 

Dial  No.  2  . 

(inches  x  10; 

Average 

Deflection  . 

(inches  x  10  ) 

75 

24 

16 

20 

135 

34 

31 

32.5 

185 

44.5 

43 

43.5 

305 

73 

72 

72.5 

400 

94 

94.5 

94 

520 

121.5 

124.5 

123 

575 

133 

139 

136 

630 

150 

151 

150.5 

760 

184 

184 

184 

805 

193 

193 

193 

865 

206.5 

205.5 

206 

1100 

264 

264 

264 

1230 

293 

292 

292.5 

l4i0 

340.5 

339-5 

340 

1505 

358 

357 

357.5 

'lot  of  data  is  shown  in  Figure  1 6 


LOAD  xIO  ,  lb 


TABLE  V.  FORCE-HEAD  SEPARATION  DATA  FOR  THE  60,  OOO-LB  BALDWIN-LIMA- 
HAMUTON  TESTING  MACHINE  MODIFIED  BY  4  LEAF  SPRINGS 


Load 

(lb  ) 

Deflection, 
Dial  No.  1  ^ 
(inches  x  10  ) 

Deflection, 

Dial  No.  2  - 
(inches  x  10  ) 

Average 

Deflection  . 

(inches  x  10  ) 

95 

27 

28 

27.5 

165 

46 

48.5 

47.2 

215 

63 

67 

65 

280 

85-5 

90.5 

88 

375 

n4 

118 

116 

470 

l4o 

145 

142.5 

VJ1 

VO 

O 

178 

185 

I81.5 

650 

193.5 

203 

197-8 

700 

211 

221 

216 

795 

237 

253.5 

245-3 

960 

290 

307 

298.5 

1070 

324 

341 

332.5 

1180 

352.5 

369 

360.8 

1270 

382 

399 

390.5 

1370 

417 

433.5 

425.3 

l48o 

454.5 

470 

462.3 

.ot  of  data  is  shown  in  Figure  17. 


LOAD 


COMPOSITE  STIFFNESS 
=  32,400  Ib/inch 


Figure  17.  Stiffness  Plot  of  60, OOO-lb  Baldvin-Lima-Hamilton 

Testing  Machine  Modified  by  Leaf  Spring  Configuration 
4l  -  Test  Series  C  and  D. 
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TABLE  VI .  FORCE-HEAD  SEPARATION  DATA  FOR  THE  60,  OOO-LB  BALDWIN-LIMA- 
HAMILTON  TESTING  MACHINE  MODIFIED  BY  3  LEAF  SPRINGS 


1  '  11  ' "  .  •  —  - 1 

Load 

Deflection 

Deflection, 

Average 

Dial  No.  1  ^ 

Dial  No.  2  . 

Deflection  ^ 

(lb  ) 

(inches  x  10  ) 

(inches  x  10  ) 

(inches  x  10  ) 

50 

18 

22 

20 

130 

44 

52 

48 

200 

70 

87.5 

78.8 

320 

115-5 

140.5 

128 

445 

162 

190 

176 

530 

193.5 

221 

207.3 

670 

247 

277 

262 

780 

290 

321.5 

305.8 

920 

343 

378 

360.5 

1045 

392 

428 

410 

13^0 

502.5 

540 

521.3 

1440 

5^3 

584 

563.5 

1540 

582 

624 

603 

Plot  of  data 

is  shown  in  Figure  38 . 

1  .  - - 

3: 


LOAD  x  10,  lb 


COMPOSITE  STIFFNESS 


Figure  18.  Stiffness  Plot  of  60,000-111  Baidu in-Lima-Hamilton 

Testing  Machine  Modified  by  Leaf  Spring  Configuration 
31  -  Tost  Series  C  and  D. 
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TABLE  VII.  FORCE-HEAD  SEPARATION  DATA  FOR  THE  60,000-LB  BALDWIN-L1MA- 
HAMUTON  TESTING  MACHINE  MODIFIED  BY  2  LEAF  SPRINGS 


Load 

(lb  ) 

Deflection, 

Dial  No.  1  r 
(inches  x  10  ) 

Deflection, 

Dial  No.  2  . 
(inches  x  10  ) 

Average 
Deflection  . 
(inches  x  104) 

45 

24 

26 

25 

105 

63 

59 

61 

210 

125.5 

129 

127.3 

290 

172.5 

167 

169.8 

380 

230.5 

229 

229.6 

470 

283.5 

286 

284.3 

585 

352 

351-5 

351.8 

695 

415 

4l4 

414.5 

780 

471.5 

471 

471.3 

920 

551-5 

552 

551.8 

1010 

604.5 

605.5 

605 

1170 

704 

706 

705 

1300 

778 

780 

779 

l44o 

861 

861 

861 

Plot  of  data  Is 

shown  in  Figure  19. 

COMPOSITE  STIFFNESS 
=16,700  It/Inch 


21  -  Tast  Series  C  and  D. 
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'ABLE  VIII.  FORCE-HEAD  SEPARATION  DATA  FOR  THE  60,  OOO-LB  BALDWIN-LIMA- 
HAMTLTON  TESTING  MACHINE  MODIFIED  BY  1  LEAF  SPRING 


Load 

(lb  ) 

Deflection 

Dial  No.  1  _ 

(inches  x  Kr) 

Deflection 

Dial  No.  2 
(inches  x  10^) 

Average 
Deflection  _ 
(inches  x  Kr) 

100 

14 

12 

13 

155 

21 

19 

20 

215 

28 

26.5 

27.3 

305 

43 

43 

43 

360 

49.5 

49.5 

49.5 

4io 

55-5 

55.5 

55.5 

480 

64.5 

64.5 

64.5 

525 

69.5 

70 

69.7 

600 

78.5 

79-5 

79 

660 

85.5 

87 

86.7 

765 

98 

100 

99 

825 

105.5 

107 

106.3 

895 

114.5 

115.5 

115 

975 

124 

125.5 

124.7 

1045 

131.5 

133.5 

132.5 

1180 

150 

151 

150.5 

1285 

162.5 

164 

163.2 

1375 

176 

177.5 

176.7 

1560 

203.5 

204.5 

204 

Plot  of  data  is  shown  in  Figure  20. 
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APPENDIX  II 

STATISTICAL  ANALYSIS  OF  THE  EXPERIMENTAL  DATA 
_ OBTAINED  IN  THE  SERIES  C  TESTS _ 

The  test  data  obtained  in  the  series  C  tests  are  presented  and  analyzed 
in  this  appendix.  Eight  levels  of  machine  stiffness  were  used  in  this 
series.  These  stiffnesses  were  obtained  as  described  in  Section  5;  the 
values  are  computed  in  Appendix  I  and  listed  in  Table  IX  together  with 
the  corresponding  buckling  loads  obtained  from  the  repeated  tests. 

A  linear  regression  analysis  is  made  on  these  data.^  For  this  purpose, 
a  regression  line  is  chosen  having  the  form 

y  =  a  +  p  (x  -  x) 

y  =  A'  +  P  (x) 


where 

y  =  buckling  load,  lb 
x  =  machine  stiffness,  lb/in. 

The  parameters  a  and  P  are  estimated  by  the  method  of  least  squares  to 
obtain  the  empirical  regression  line 

Y  =  a  +  b(x-x)  =  A+bx 


From  the  data  presented  in  Table  X,  the  mean  buckling  load  may  be  computed 


as 


456.72  lbs. 
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and  for  this  case  is  given  as 


572057 

259,310,973,000 


-  0.2206  x  10 


5 


in. 


Thus,  the  enquiries!  regression  line  becomes 

Y  =  456.72  -  0.2206  x  10"5  (x  -  x) 

This  line  is  drawn  In  Figure  12  through  the  discrete  data  points, 

Since  more  them  one  value  of  buckling  load  was  determined  at  each  stiff¬ 
ness,  the  hypothesis  regarding  the  linearity  of  the  regression  curve  may 
be  tested  by  comparing  the  variation  about  the  regression  line  with  that 
existing  within  sets.  Linearity  is  rejected  if 

2 

sy 

F  =  —  a  Fa;  n-2,  n(k  -  l) 

S1 

where  n  =  number  of  sets 

k  =  number  of  values  within  each  set 


p 

-Yi)/(n-2)  (variance  about 
regression  line) 


s 


2 

1 


I 


71=1 


„  2 

(y.-n  -  y. )  /  n(k  -  l)  (variance 

within  sets) 


Fq;  n-2.  n(k  -  l)  =  the  100  Of  percentage  point  of  the  F  distribution 
with  (n-2)  and  n(k-l)  degrees  of  freedom 

y^  =  the  T)  th  value  in  the  i  th  set 

yi  =  the  average  value  of  the  i  th  set 

Y  =  the  empirical  value  of  y  corresponding  to  x 
n  2 

It  is  evident  that  ^  (y^  -  Y^)  is  Just  the  sum  of  squares  of 

i  =  1 

the  deviations  of  the  average  value  of  each  set  about  the  fitted  line. 
The  computations  are  carried  out  in  Table  X  using  the  relation 
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1  (yi  ■  Yi)2  =  L  (yi  ■ y)2 


i  =  i 


i  =  i 


1  2 


£  (x±  -  x)  (y±  -  y) 


i  =  1 


I  (*i  =  x)2 


i  =  1 


T}jus,  variance  of  y  about  the  mean  a  +  b 
a  ,  is  estimated  from  the  data  as 

Sg  =  6.28  - 


(x  -  x  ),  described  exactly  by 

(572,057 )2  2 

259,310,973, ooo  =  5-02  lb 


Likewise,  it  is  evident  that  «  —  \2  .  ..  _  , 

v  y  (y^  -  y±)  is  the  sum  for  all  sets 

i=l  T|=l 


of  the  sum  of  the  squares  of  the  deviations  cf  the  individual  values 
within  each  set  about  the  mean  value  of  the  set.  The  computations  are 
carried  out  in  Table  X  using  the  relation 


l 


i=l 


I  Wti  -  -  I 

71=1  i=l 


l 


71=1 


K/Ai 


71=1 


2 

Before  computing  s^,  the  homogeneity  of  the  several  set  variances  is 
tested  by  means  or  Cochran's  test.  The  hypothesis  of  equality  is 
accepted  if 


where 


g  ■  S1  *  tg 

n 

i=l 

k 

s?=  I  ^yi71  "  yi^  /  ^ki  “  ^ 

1  7]  =  1 
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is  the  ith  fl£t  variance  and  depends  upon  the  level  of  significance 
a,  n  and  k.  Choosing  a  =  5  percent,  is  O.36O.  Then,  from 

Table  IX,  \  *  J 

g  =  ^  =  0.202  <  O.36O 

Thus,  the  eight  variances  do  not  differ  significantly  and  nay  be  combined 
to  give  the  single  estimate 


s 


2  _  218 
1  “  To 


=  5.^5  lb.2 


It  is  found  that 

FQ;  n-2,  n(k-l)  =  F0<05j6,4O  =  2.34 

Substitution  yields 

F  *  ini  -  °'921 

Thus,  the  hypothesis  of  linearity  is  accepted  at  the  95 “percent  level. 
This  leads  to  a  pooled  estimate  of  variance  about  the  regression  line: 

i  -  V2  +  1  1  <ytn  ■ 

s  =  1=1 _ i=l  T]=l _ 

nk  -  2 

s2  .  6 (;-°2)  ♦  21a  =  5.3939 

and  a  standard  deviation: 


s  =  2.3225 


The  variance  of  the  estimate  b  of  the  slope  and  the  estimate  A  of  the 
intercept  sure  normally  distributed  with  means  6  and  A',  respectively.  The 
variance  of  b  is  given  by 


From  Table  X  end  the  estimate  for 


2 

0  > 


we 


may  write 
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and  for  A'  by 

A  -  ^x^jkn  -  2  A 

where  2  iS  the  °/2  percentage  point  of  Student's  t 

distribution.  Choosing  a  =  5  percent,  the  95-percent  confidence  interval 
estimate  for  8  is 

b  +  t  025.i^CTb  =  [-0.2206  x  10"5  +  2.015(0.1662)  x  10"5] 
=  [-0.5958;  0.1546]  x  10"5  in. 

The  confidence  interval  for  A'  is 


4l 


A  -  t.025;  h6  °A  =  [  456,97  -  2-°15(0.6996)  ] 
=  [  455.56;  458.38  ]  lb 


The  theoretical  slope  £  and  the  intercerpt  A*  are  examined  statistically 
to  check  the  initial  assumption  that  buckling  load  is  a  linear  function 
of  machine  rigidity.  The  significance  tests  for  these  coefficients  are 
given  in  Table  XI. 


The  empirical  regression  line  shown  in  Figure  12  is  almost  horizontal. 

In  addition,  the  95-percent  confidence  intervc1  for  P  includes  the  possib¬ 
ility  of  a  zero  slope.  This  suggests  that  there  is  no  relationship 
between  x  (machine  stiffness)  and  the  mean  value  of  y  (mean  buckling  lord), 
and  that  the  small  empirical  slope  b  is  due  to  accidental  variation  of  the 
data . 


This  hypothesis  is  tested  by  putting  P  =0  in  Table  XI  .  The  test 
statistic  is 


t 


0.2206 

0.1862 


-1.185 


At  the  5-percent  level  of  significance,  the  criterion  for  rejection 
becomes 

I  1 1  21  ^ytejkn  -  2  =  t0.025;46  =  2’015 

Hence 

1.185  <  2.015 


The  hypothesis  that  the  theoretical  slope  is  zero  is  accepted  at  the 
95 -percent  level. 

A  measure  of  the  sensitivity  of  the  analysis  is  obtained  from  an  examina¬ 
tion  of  the  0C  curve  at  this  level.  For  a  sample  size  of  48,  there  is  a 
95-percent  probability  of  detecting  a  value  of  d  =  0.6. 


For  this  test. 


0.6 


Hence, 

|  |  =  0.6  (  0.1862  x  10"5)  ^”47" 

|  P]_  |  =  0.7659  x  10~5 
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or 


Thus,  the  hypothesis  -p  =  0  would  be  rejected  at  the,-95 -percent  level 
if  it  differed  from  zero  by  as  little  as  O.7659  x  10"')  in.  A  slope  of 
this  magnitude  may  be  seen  in  better  perspective  if  we  note  that  it  implies 
a  change  at  zero  stiffness  given  by 

\  -  V 

=  (0.7659  X  io'5)(m,275) 

=  0.85  lb 

which  is  only 

(100)  O.83  =  0.23$ 

J+56.72 


of  the  average  critical  load. 

The  assumption  that  no  relationship  exists  between  x  and  the  mean  value 
of  y  may  be  tested  further  by  assuming  that  the  theoretical  intercept  A* 
is  equal  to  the  mean  value  of  the  buckling  load.  The  95-percent  confi¬ 
dence  interval  for  A1  includes  the  possibility  of  a  value  equal  to  y. 

To  test  this  hypothesis,  A1  is  set  equal  to  y=456.72  in  Table  XX.  The 
test  statistic  is  0 


Since 


t  =  A  -  A1 

0 


[456.72  +  (0.2206  x  10"^)  (111,275)  -  456.72 
t  =  U  0.3940 


t 


0.2455 

0.3940 


0.6231 


O.623  <  t 


.025; 46 


2.015 


the  hypothesis  -A’  =  y  is  accepted  at  the  95-percent  level. 


Then,  from  the  0C  curve. 


Hence, 


d  =  ,j  A'q  -  A'j  =  0.6 
CA  lj  kn-1 

I  A'o  '  A'i  i  =  0*6(0.3940) 

I  A*  -  A'.  I  =  1.62  lb 
1  o  1  1 
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Thus,  the  hypothesis  -A'  =  y  would  be  rejected  at  the  95_percent  level 
if  the  theoretical  intercept  differed  from  the  mean  value  by  as  little  as 
1.62  lb,  or  only  [lOO(l£2)/456.72]  =  0.35$  of  the  average  critical  load. 


TABLE  IX.  LINEAR  REGRESSION  ANALYSIS  -  SERIES  C  RESULTS 
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TABLE  X.  LINEAR  REGRESSION  ANALYSIS  -  SERIES  C  TESTS 
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TABLE  XI. 

SIGNIFICANCE  TEST  FOR  SLOPE  AND  INTERCEPT  OF  STRAIGHT  LINE 
(See  Reference  11 ) 

Hypothesis  Test 

Statistic 

Criteria  for 
Rejection 

Operating 
Characteristic 
Abscissa  Value 

t 

1  t  |  * 

d 

e  =  „o 

b-PQ 

V^jkn  -  2 

1  -  h  1 

CTb 

°b  /kn  -  1 

> 

II 

> 

A  -  A' 

o  0 

*0/2; kn  -  2 

1  A'  -  A\  1 

1  o  1  1 

CTA  /kn  -  1 

n 

=  number  of  sets 

k 

=  number  of  values  within  each  set 

a 

=  level  of  significance 

j  P  and  A' 
o  o 

=  hypothetical  values 

of  slope  and  intercept,  respectively 

P.^  and  A'^ 

i 

! 

=  the  variations  of  slope  and  intercept  from  the  values 
which  can  be  detected  for  the  given  sample  size  and 
chosen  probability.  In  finding  the  value  of  d,  the 
(n  -  l)  curve  in  Reference  11  is  used. 

APPENDIX  III 

PRESENTATION  AND  STATISTICAL  ANALYSIS 
OF  THE  EXPERIMENTAL  DATA  OBTAINED 
_ IN  THE  SERIES  D  TESTS 


The  test  data  obtained  in  the  series  D  tests  are  presented  and  analyzed 
in  this  appendix.  Eight  levels  of  machine  stiffness  were  used  in  this 
series.  These  stiffnesses  were  obtained  as  described  on  page  3  j  the 
values  are  computed  in  Appendix  I  and  are  listed  in  Table  XII  together 
with  the  corresponding  buckling  loads  obtained  from  the  repeated  tests. 

A  linear  regression  analysis  is  made  on  these  data.^  For  this  purpose, 
a  regression  line  is  chosen  having  the  form 

y  =  Of  +  P  (x  -  x) 

y  =  A'  +  p  (x) 

where  y  =  buckling  load,  lb 

x  =  machine  stiffness,  lb  /in. 

The  parameters  a  and  (3  are  estimated  by  the  method  of  least  squares  to 
obtain  the  empirical  regression  line 

Y  =  a  +  b  (x  -  x)  -  A  +  bx 


From  the  data  presented  in  Table  XIII, 
computed  as 


a 


the  mean  buckling  load  may  be 


=  ?00.7  =  112.59  lb 
£3 


while  the  mean  value  of  machine  stiffness  used  is  given  by 

n 

x  =  i=i  =  890200  =  111,275  lb  / in. 

— 8 —  — 8” 

The  slope  b  is  defined  by 

n 

k  I  (xi  ■  x)  ■  y) 

b  =  i=l _ 

(xi  "  x)2 

i=i 


and  for  this  case  is  given  as 


611, jkk 

259,310,973,000 

kQ 


-0.23576  x  10"5 


in. 


Thus,  the  empirical  regression  line  becomes 


-5 

Y  =  112.59  -  0.23576  x  10  (x  -  x) 

This  line  is  drawn  in  Figure  13  through  the  discrete  data  points,  y^. 

Since  more  them  one  value  of  buckling  load  was  determined  at  each  stiff¬ 
ness,  the  hypothesis  regarding  the  linearity  of  the  regression  curve 
may  be  tested  by  comparing  the  variation  about  the  regression  line  with 
that  existing  within  sets.  Linearity  is  rejected  if 

2 

F  =  >  Fa;  n-2,  n(k-l) 

81 

where 

n  =  number  of  sets 


k  -  number  of  values  within  each  set 


s^  =  k 


2  =  k  \  (y .  -  Y, )  /  (n  -  2)  (variance  about 

L  1  1  regression  line) 


i=l 

n  k 


31  =  £  ^  (yiTj  -  y^2/  n(k  -  l)  (variance  within 

11.1  sets) 


TJ» 

a;  n  -  2,  n(k  -  l)  =  the  100  a  percentage  point  of  the  F 

distribution  with  (n  -  2)  and  n(k  -  l) 
degrees  of  freedom 

yiTj  =  the  7)th  value  of  the  i  th  set 

y^  =  the  average  value  of  the  i  th  set 

Y.,  =  the  empirical  value  of  y  corresponding  to  x. 

j>  1 

It  is  evident  that 

i  -  v2 

i=l 

is  Just  the  sum  of  squares  of  the  deviations  of  the  average  value  of 
each  set  about  the  fitted  line.  The  computations  are  carried  out  in 
Table  XIII  using  the  relation 
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1  <yi-V2  ■  l  <*i-y)2- 


i  =  1 


i  =  1 


n 


£  (xi  ■  x)(yt  -  y) 

Li  =  1 


I  (xi  -  x' 


i  =  1 

Thus,  the  variance  of  y  about  the  mean  OL  +  0(x  -  x),  described  exactly 
by  o,  is  estimated  from  the  data  as 


s2  =  7.89 _ 

b2  259,310,973,000 


=  6.45  lb 


Likewise,  it  is  evident  that 


I  1  (ym  ’  y±)2 


i=l  Tl=l 


ill 


is  the  sum  for  all  sets  of  the  sum  of  the  squares  of  the  deviations  of 
the  individual  values  within  each  set  about  the  mean  value  of  the  set. 
The  computations  are  carried  out  in  Ttble  XII  uBing  the  relation 

i  £  (yiri  - y/  ■  £  i  yi2  -  £  i  <yi//ki 

i=l  T|=l  i=l  71=1  i=l  71=1 

2 

Before  computing  a. ,  the  homogene i|^  of  the  several  set  variances  is 
tested  by  means  oxCochran 1  s  test.  The  hypothesis  of  equality  is 
accepted  if 

g  =  lareeBt  °i  s  * 
n 

H 

i=l 
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where 


k 


1  (ym  -  yi)2/  ^ki 
Thl 


1) 


is  the  i  th^£et  variance  and  depends  upon  the  level  of  significance 
Of,  n  and  k.  Choosing  a  =  5  percent,  z  o\  is  O.36O.  Then,  from 
Table  XII, 


g  = 


6.2 

29.2 


=  0.212  <  0.360 


Thus,  the  eight  variances  do  not  differ  significantly  and  may  be  combined 
to  give  the  single  estimate 


2  146 

8l  =  To  = 


3.65 


It  is  found  that 


Fa;n-2,  n(k-l)  =  F0.05;6,60  =  2‘-^ 

Substitution  yields 

F  =  =  I.76  <  2.34 


Thus,  the  hypothesis  of  linearity  is  accepted  at  the  95 -percent  level. 
This  leads  to  a  pooled  estimate  of  variance  about  the  regression  line : 


2 

s 


I  frm  -  *i>2 

M _ : _ 


2  .  6(6.^  +  M  ,  4.0152 


and  a  standard  deviation: 


s  =  2.004  lb 


The  variance  of  the  estimate  b  of  the  slope  and  the  estimate  A  of  the 
intercept  are  normally  distributed  with  means  0  and  A1,  respectively. 
The  variance  of  b  is  given  by 


2 

_ o _ 

n 

kI  <*!  ‘  k>2 

i=l 
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-  ,*■**<» 


From  Table  XIII  and  the  estimate  for  a  ,  we  may  write 


CTb  “  6(259,310,973,000)  =  2-58o605  x  10 


Thus,  the  standard  deviation  is  given  as 


ab  =  0.1606  x  10 


-5 


The  variance  of  A  is  given  by 

2  2 
aA=  a 


J.  + 

nk 


-2 

x 


n 
\ 

k 

i 

i=l 


(xi  -  x)‘ 


and  by  substituting  appropriate  values  from  Table  XIII  we  arrive  at 

=  0.115604 


Then, 


- 

^^.0152  | 


12,382,125,600 
TE  +  6(259,310,913,000) 


aA=  0.3400 


In  addition  to  the  point  estimates  of  slope  and  intercept,  confidence 
intervals  for  p  and  A'  may  be  established  with  confidence  coefficient 
1  -  a.  They  are  given  for  P  by 


”  ±  Vs  -  2'’'’ 


and  for  A'  by 

A  -  ta/2;  kn  -  2aA 

where  t a^.  g  Is  the  100  ce/2  percentage  point  of  Student's  t 

distribution.  Choosing  a  =  5  percent,  the  95 -percent  confidence  interval 
estimate  for  P  is 

b  ±  t.025j  46ab  =  [-°-2358  x  10-5 

=  [-0.559;  0.088] 

The  confidence  interval  for  A'  is 


+  2.015(0.1606)  x  10-5] 
x  10" 5  in. 
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A  ±  t  025.  46oa=  [  112.85  +  2.015(0.3^*00)] 
=  [  112.17;  113.5;’  ]  lb 


The  theoretical  slope  p  and  the  intercept  A'  are  examined  statistically 
to  check  the  initial  assumption  that  buckling  load  is  a  function  of 
machine  rigidity.  The  significance  tests  for  these  coefficients  are 
given  in  Table  XI. 


The  empirical  regression  line  shown  in  Figure  13  is  almost  horizontal. 

In  addition,  the  95 -percent  confidence  interval  for  p  includes  the 
possibility  of  a  zero  slope.  This  suggests  that  there  is  no  relationship 
between  x  (machine  stiffness)  and  the  mean  value  of  y  (mean  buckling 
load),  and  that  the  small  empirical  slope  b  is  due  to  accidental 
variation  of  the  data. 


This  hypothesis  is  tested  by  putting  P  =  0  in  Table  XI.  The  test 
statistic  is 


t 


_b  _  _  0.23476 

a  ~  0.1606 

b 


-1.468 


At  the  5-percent  level  of  significance,  the  criterion  for  rejection 
becomes 


I  t  I  2  fysikn  -  2  =  t0.025;  46  =  2,015 


Since 


1.468  <  2.015 


the  hypothesis  that  the  theoretical  slope  is  zero  is  accepted  at  the  95- 
percent  level. 

A  measure  of  the  sensitivity  of  the  analysis  is  obtained  from  an  exam¬ 
ination  of  the  OC  curve  at  this  level.  For  a  sample  size  of  48,  there 
is  a  95-percent  probability  of  detecting  a  value  of  d  =  0.6. 


For  this  test, 

d  =  ^o  "  =  0.6 

°b  /kn  -  1 

Hence, 

|  P1  |  =  0.6  (0.1606  x  10"5)  AT 

|  Px  j  =  0.6606  x  10'5 
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Thus,  the  hypothesis  -0  =  0  would  be  rejected  at^the  95-percent  level  if 
it  differed  from  zero  by  as  little  as  0.6606  x  lO-'*  in.  A  slope  of  this 
magnitude  may  be  seen  in  better  perspective  if  we  note  that  it  implies  a 
change  at  zero  stiffness  given  by 


Ay  =  Pxx 

=  (0.6606  x  10  5)  (111275) 
=  0.74  lb 

which  is  only 

(100)  0.74  =  0.66# 

11.259 


of  the  average  critical  load. 

The  assumption  that  no  relationship  exists  between  x  and  the  mean  value  of 
y  may  be  tested  further  by  assuming  that  the  theoretical  intercept  A' 
is  equal  to  the  mean  value  of  the  buckling  load.  The  95 -percent  confidence 
interval  for  A'  includes  the  possibility  of  a  value  equal  to  y.  To  test 
this  hypothesis,  A'  is  set  equal  to  y  =  112.59  in  Table  XI.  The  test 
statistic  is 


t  =  A  -  A’o 


t  _  [112.59  +  (0.23576  X  10"5)  (111275) 

0.3^00 


-  112.59 


,9] 


0.2623  _  n 
t  =  o735oo  =  °’(n 5 


Since 


0.7715  «  t  =  2.015 

the  hypothesis  -A'  -  y  is  accepted  at  the  95-percent  3,evel. 
Then,  from  the  OC  curve, 


Hence, 


aA  /kn  -  1 

I  A'  -  A’  |  =  0.6(0.3400)  JW 
1  o  1 

I  A'  -  A*  I  =  1.40  lb 
1  o  1  1 
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f  .n*  - 


Thu 3,  the  hypothesis  -A'  =  y  would  he  rejected  at  the  95-Percen^  lsvel 
if  the  theoretical  intercept  differed  from  the  mean  value  by  as  little  as 
1.40  lb  or  only 

io°(i>£L . 

112.59 

of  the  average  critical  load. 


TABLE  XII.  LINEAR  REGRESSION  ANALYSIS  -  SERIES  D  RESULTS 
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TABLE  XIII.  LINEAR  REGRESSION  ANALYSIS  -  SERIES  D  RESULTS 
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.  Figure  5(c) 


